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We report on a search for gravitational waves from the coalescence of compact binaries during the third 
and fourth LIGO science runs. The search focused on gravitational waves generated during the inspiral phase 
of the binary evolution. In our analysis, we considered three categories of compact binary systems, ordered by 
mass: (i) primordial black hole binaries with masses in the range 0.35Mq < mi, rri2 < I.OMq, (ii) binary 
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neutron stars with masses in the range I.OMq < mi,m2 < 3.0Mq, and (iii) binary black holes with masses in 
the range 3.0M© < mi, 7712 < m ma x with the additional constraint mi + mi < m max , where m max was set 
to 4O.OA/0 and 8O.OA/0 in the third and fourth science runs, respectively. Although the detectors could probe 
to distances as far as tens of Mpc, no gravitational-wave signals were identified in the 1364 hours of data we 
analyzed. Assuming a binary population with a Gaussian distribution around O.75-O.75M0, 1.4-1.4Mq, and 
5.O-5.OM0, we derived 90%-confidence upper limit rates of 4.9 yr _1 L^" 1 for primordial black hole binaries, 
1.2 yr _1 L]~ for binary neutron stars, and 0.5 yr _1 L]~ for stellar mass binary black holes, where L10 is 10 10 
times the blue light luminosity of the Sun. 

PACS numbers: 95.85.Sz, 04.80.Nn, 07.05.Kf, 97.80.-d 



I. OVERVIEW 

While gravitational radiation has not yet been directly de- 
tected, observations of the orbital decay of the first binary pul- 
sar PSR B1913+16 fl 0] have provided significant indirect 
evidence for their existence since the late eighties. Indeed, 
observations have revealed a gradual inspiral to within about 
0.2 percent of the rate expected from the emission of gravita- 
tional radiation |3J]. As orbital energy and angular momentum 
are carried away by gravitational radiation, the two compact 
objects in a binary system become more tightly bound and or- 
bit faster until they eventually merge. The gravitational wave 
signals emitted by the merging of binary systems made of pri- 
mordial black holes, neutron stars, and/or stellar mass black 
holes can be detected by ground-based detectors. The detec- 
tion rate depends on the merger rate, which in turn depends 
on the rate of ongoing star formation within LIGO's detection 
volume, described in greater detail in ||4j] and as measured 
by the net blue luminosity encompassed in that volume (see 
Sec.|IVj. 

Several direct and indirect methods can be applied to infer 
the mergerrate expected per unit Liq, where Liq is 10 10 times 
the blue solar luminosity. Merger rates for binary neutron 
star (BNS) systems can be directly inferred from the four sys- 
tems observed as binary pulsars that will merge in less than a 
Hubble time; the basic methodology was originally applied by 
10 0]. The current estimates based on all known BNS suggest 
that the merger rate lies in the range 10-170 x 10 _6 yr _1 L 10 1 
10,0]. This range is at 95% confidence for a specific model of 
the Galactic population (model #6 in the references), which 
represents our current understanding of the radio pulsar lu- 
minosity function and their Galactic spatial distribution. The 
most likely rate for the same model is 50 x 10- 6 yr _1 £io 
|0 0]. The estimated BNS merger rate makes the detection 
of a signal from such an event unlikely, though possible, with 
the current generation of gravitational-wave detectors. In con- 
trast, there is no direct astrophysical evidence for the exis- 
tence of binary black hole (BBH) or black hole/neutron star 
binaries, but they are predicted to exist on the basis of our 
current understanding of compact object formation and evo- 
lution. The search for gravitational waves emitted by BBH 
systems is particularly interesting since it would provide di- 
rect observation of these systems. Merger rate estimates are 
currently obtained from theoretical population studies of bi- 
naries in galactic fieldsJ^JT^JTll [H [H 0, [ij [H] or in 
dense stellar clusters Il7l 18T 19ll. Because these studies dif- 
fer significantly in their assumptions and methodology, it is 



difficult to assess all the literature and assign relative likeli- 
hoods to merger different merger rates for black hole binaries. 
However, in the case of field binaries, estimates for the rela- 
tive likelihood can be obtained by widely exploring several of 
the parameters of the population models, while ensuring those 
models reproduce the BNS merger rates derived from the ob- 
served sample 1I20I l2lll . Based on this study, the merger rates 
for BBH and black hole/neutron star binaries are found to lie 
in the ranges (at 95% confidence) 0.1 — 15 x 10~ 6 yr _1 L^ 1 
and 0.15 — 10 x 10~ 6 yr -1 !^ respectively, with most likely 
merger rates of 0.6 x 10 -6 yr and 1.3 x 10~ 6 yr _1 L^ L0 1 . 

Although drawn from a single study, the simulations cover 
such a uniquely wide parameter space that these rate ranges 
are consistent with the existing literature on BBH and black 
hole/neutron star merger rates. It has also been discussed 
in the literature that some fraction of all dense clusters may 
form many inspiraling BBH; although the current rate pre- 
dictions are considered highly uncertain and the systematic 
uncertainties are not yet understood, rates as high as a few 
events per year detectable by initial LIGO have been reported 
lfl7l fla, [Hal . Furthermore, indirect evidence suggests that 
short, hard gamma-ray bursts (GRB)s could be associated 
with the coalescence of a BNS or a black hole/neutron star 
binary. Recent estimates suggest that the rates of these events 
could be in excess of about 1 x 10~ 6 yr -1 /^ 1 |22tl . There 
may also exist sub-solar-mass black hole binary systems, with 
component objects that could have formed in the early uni- 
verse and which contribute to galactic dark matter halos l23ll ; 
we refer to such lower-mass compact binary coalescences as 
primordial black hole (PBH) binaries. 

The Laser Interferometer Gravitational-wave Observatory 
(LIGO) Scientific Collaboration (LSC) operates four interfer- 
ometric detectors. Three of these are from the U.S. LIGO 
project j24ll25Tl . two of them, with 4 km and 2 km long arms, 
are co-located in Hanford, WA (called HI and H2, respec- 
tively) and a third detector, with 4 km long arms, is located 
in Livingston, LA (called LI). The LSC also operates the 
British-German GEO 600 detector j26Tl . with 600 m long arms 
that is located near Hannover, Germany. Only data from the 
LIGO detectors were used in this analysis, however, due to the 
relative sensitivity of the detectors. 

We report on a search for gravitational waves emitted by 
coalescing compact binaries in the data taken by the LIGO 
detectors in late 2003 (Oct 31, 2003-Jan 9, 2004) and early 
2005 (Feb 22, 2005-March 24, 2005) which correspond to the 
third (S3) and fourth (S4) science runs, respectively. During 
S3 and S4, the LIGO detectors were significantly more sensi- 
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tive than in our previous science runs 1T271 l28i I29I l30ll . This 
improvement can be quantified in terms of the inspiral horizon 
distance of each detector which is defined as the distance at 
which an optimally located and oriented binary system would 
give expected signal-to-noise ratio (SNR) equal to 8. For in- 
stance, HI, the most sensitive detector during S4, had hori- 
zon distance averaged over the duration of the run of 5.7 Mpc, 
16.1 Mpc, and 77.0 Mpc, for a 0.5-0. 5A/ Q , 1.4-1.4M , and 
lO-lOM© systems, respectively. Consequently, during S3 and 
S4, the detectors were sensitive enough to detect inspiral sig- 
nals from hundreds of galaxies as shown in Fig.Q] 

The paper organization is as follows. In Sec. [II] we briefly 
describe the data analysis pipeline and present the parameters 
used in the S3 and S4 science runs. In particular, Sec. Ill CI 
describes the division of the search into 3 categories of bina- 
ries: PBH binary, BNS, and BBH inspirals. In Sec. [ml we 
present the results of the search, including the accidental rate 
estimates and loudest candidates found from the different sci- 
ence runs and categories of binary systems that we considered. 
Finally, Sec.lIVIdescribes the upper limits set by this analysis. 

H. THE DATA ANALYSIS PIPELINE 

The analysis pipeline used to search the S3 and S4 data re- 
ceived substantial improvements over the one used in our pre- 
vious searches J27il28il29l[30ll . The pipeline is fully described 
in a set of companion papers lf3ll [52, HH; this section intro- 
duces the aspects of our analysis methodology that are essen- 
tial to comprehend the search and the final upper limit results. 
We emphasize the differences between the BBH search and 
the PBH binary /BNS searches. 



A. Coincident data and time analyzed 

The first step of the analysis pipeline is to prepare a list 
of time intervals represented by a start and end time, during 
which at least two detectors are operating nominally. Requir- 
ing coincident signals from two or more detectors reduces the 
accidental rate by several orders of magnitude and increases 
our detection confidence. 

In S3, we required both Hanford detectors to be operating; 
analyzed times belonged either to triple H1-H2-L1 or double 
H1-H2 coincident times. In S4, times when HI was operating 
but H2 was not (and vice-versa) were also analyzed, there- 
fore all permutations of double coincident times were pos- 
sible, in addition to the triple coincident times. The break- 
down of times analyzed, common to all searches, is given in 
Table [I] A fraction of these times (about 9%), playground 
times, was used to tune the search parameters. This tuning 
was performed in order to suppress background triggers orig- 
inating from instrumental noise so as to efficiently detect the 
gravitational wave signals (measured using simulated injec- 
tions, as described in Section IIII Ab . In order to avoid po- 
tential bias, upper limits (Section |IV| ) are derived using the 
non-playground data only. However, candidate detections are 
drawn from the full data set. 




10 Total binary mass (Mq) 



FIG. 1: Blue-light luminosities and horizon distances for LIGO's 
observatories. In the left panel, the horizontal bars represent the non- 
cumulative intrinsic blue-light luminosity of the galaxies or clusters 
within each bin ofphysical distance, as obtained from a standard as- 
tronomy catalog U. Some bins are identified by the dominant con- 
tributor galaxy or cluster. The merger rate of binaries within a galaxy 
or cluster is assumed to scale with its blue-light luminosity. The de- 
tectability of a binary depends on the effective distance between the 
source and detector, which is dependent on both the physical dis- 
tance separating them and their relative orientation (see Eq.[2]l. The 
solid line shows the cumulative blue-light luminosity as a function of 
effective distance (hereafter, the effective cumulative blue-light lumi- 
nosity), of the binary sources which would be observed by the LIGO 
detectors if they had perfect detection efficiency (i.e., all binaries are 
detectable). Explicitly, a binary in M31 (at a physical distance of 
0.7 Mpc) with an effective distance of 5 Mpc will contribute to the 
intrinsic luminosity at 0.7 Mpc and will contribute to the effective 
cumulative luminosity at a distance of 5 Mpc. Although a binary 
will have slightly different orientation with respect to each LIGO 
observatory and therefore slightly different effective distances, the 
difference in the effective cumulative luminosity shown on this plot 
would not be distinguishable. The effective cumulative luminosity 
starts at 1.7 L10 (Milky Way contribution), and begins increasing 
at a distance of ~ 1 Mpc, with the contribution of nearby galax- 
ies M31 and M33. The cumulative luminosity observable by our 
search (not shown), as expressed in Eq. [8] depends also on the de- 
tection efficiency of our search (see Fig.[5j and will be less than the 
effective cumulative luminosity. In the right panel, the curves rep- 
resent the horizon distance in each LIGO detector as a function of 
total mass of the binary system, during S3 (dashed lines) and S4 
(solid lines). We also plot the horizon distance of LI during S2. The 
sharp drop of horizon distance around a total mass of 2Mq is related 
to a different lower cut-off frequency, /l, used in the PBH binary 
search and the BNS/BBH searches. The /_l values are summarized 
in Table Inl The high cut-off frequency occurs at the last stable orbit. 
The horizon distance for non equal mass systems scales by a factor 
(mi + 7712). 



We compiled a list of time intervals when the detectors had 
poor data quality 1I32I l34ll . In S3, this selection discarded 5% 
of H1/H2 as a result of high seismic noise and 1% of LI 
data as a result of data acquisition overflow. In S4, 10% of 
H1/H2 data was discarded mostly due to transients produced 
when one Hanford detector was operating but the other was 
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TABLE I: Times analyzed when at least two detectors were oper- 
ating. The times in parentheses exclude playground times, which 
represents about 9% of the data and is used to tune the search. 





S3 


S4 


H1-H2-L1 times 


184 (167) hrs 


365 (331) hrs 


H1-H2 times 


604 (548) hrs 


126 (114) hrs 


HI -LI times 




46 (41) hrs 


H2-L1 times 




39 (35) hrs 


Total times 


788 (715) hrs 


576 (521) hrs 



not. A gravitational wave arriving during one of the vetoed 
times could, under certain conditions, still be detected and val- 
idated. However, neither playground times nor vetoed times 
are included when computing the upper limits presented in 
Section HV] 

B. Filtering 

In the adiabatic regime of binary inspiral, gravitational 
wave radiation is modeled accurately. We make use of a va- 
rietyof approximation techniques l35l [36j [37j [38j [39j |40j |4li 
I42L I43I1 which rely, to some extent, on the slow motion of the 
compact objects which make up the binary. We can represent 
the known waveform by 

h(t) = ^^A(t) cos Wt)~M (1) 

where cp is some unknown phase, and the functions A(t) and 
cf)(t) depend on the masses and spins of the binary. Although 
spin effects can be taken into account l44ll . they are estimated 
to be negligible over much of the mass range explored in this 
search and will be neglected here. Since the gravitational 
wave signal we are searching for is known, the matched fil- 
tering method of detection constitutes the cornerstone of our 
analysis. In both PBH binary and BNS searches, we use 
physical template families based on second order restricted 
post-Newtonian waveforms in the stationary-phase approxi- 
mation ll37l, l45ll . In the BBH search, we use a phenomenolog- 
ical template family ll46tl so as to palliate uncertainties in the 
gravitational-wave templates, which become significant in the 
LIGO band for higher mass systems. The template matched 
filtering will identify the masses and coalescence time of the 
binary but not its physical distance D. The signal ampli- 
tude received by the detector depends on the detector response 
functions F + and F x , and the inclination angle of the source t, 
which are unknown. We can only obtain t he ef fective distance 
D e g, which appears in Eq. ([TJ defined as II47I1 : 

D cS = = = • (2) 

JFI(1 + cos 2 i) 2 /4 + Fl (cos t) 2 

The effective distance of a binary may be larger than its phys- 
ical distance. 



C. Inspiral search parameters 

We searched for PBH binaries with component masses be- 
tween 0.35 Mq and 1 Mq, and BNS with component masses 
between 1 Mq and 3 Mq. We also searched for BBH systems 
with component masses between 3 Mq and m max , where 
m max was set to 40 Mq and 80 Mq in S3 and S4, respec- 
tively. In addition, the total mass of the systems was also 
constrained to be less than m max . The larger mass range 
in S4 is due to improvement of the detector sensitivities at 
low frequency. This classification of binaries into three cat- 
egories was driven primarily by technical issues in the data 
analysis methods. In particular, the waveforms differ signif- 
icantly from one end of the mass scale to the other: gravita- 
tional waves from lower mass binaries last tens of seconds in 
the LIGO band and require more templates to search for them, 
as compared to the higher mass binaries (see Table Hill. 

For each search, we filtered the data through template 
banks designed to cover the corresponding range of compo- 
nent masses. The template banks are generated for each de- 
tector and each 2048-second data stretch so as to take into 
account fluctuations of the power spectral densities. In the 
PBH binary and BNS searches, the algorithm devoted to the 
template bank placeme nt [l33ll is identical to the one used in 
previous searches 11281 12911 . In the BBH search, we used a 
phenomenological bank placement similar to the one used in 
the S2 BBH search [f30h - The spacing between templates gives 
at most 5% loss of SNR in the PBH binary and BBH banks, 
and 3% in the BNS bank. The average number of templates 
needed to cover the parameter space of each binary search are 
shown in Table [II] and are indicative of the relative computa- 
tional cost of each search. 

TABLE II: The target sources of the search. The second and third 
columns show the mass ranges of the binary systems considered. The 
fourth column provides the lower cut-off frequency, /l, which set 
the length of the templates, and the fifth column gives the average 
number of templates needed, iVb. The last column gives the longest 
waveform duration, T max . 





m min (M Q ) 


m max (M ) 


/l(Hz) 


N b 


2~max(s) 


S3,S4 PBH 


0.35 


1.0 


100 


4500 


22.1 


S3 BNS 


1.0 


3.0 


70 


2000 


10.0 


S4 BNS 


1.0 


3.0 


40 


3500 


44.4 


S3 BBH 


3.0 


40.0 


70 


600 


1.6 


S4 BBH 


3.0 


80.0 


50 


1200 


3.9 



For each detector, we construct a template bank which we 
use to filter the data from the gravitational wave channel. Each 
template produces an SNR time series, p{t). We only keep 
stretches of p(t) that exceed a preset threshold (6.5 in the PBH 
binary and BNS searches and 6 in the BBH case). Data reduc- 
tion is necessary to cope with the large rate of triggers that 
are mostly due to noise transients. First, each SNR time se- 
ries is clustered using a sliding window of 16 s as explained in 
ll3lll . Then, surviving triggers from all templates in the bank 
are clustered, so that only the loudest template trigger is kept 
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in fixed intervals of 10 ms (PBH binary and BNS) or 20 ms 
(BBH). These triggers constitute the output of the first inspi- 
ral filtering step. To further suppress false triggers, we require 
additional checks such as coincidence in time in at least two 
detectors, as described below. 



D. Coincidence parameters and combined SNR 

In the PBH binary and BNS searches, we require coinci- 
dence in time, chirp mass A4 C = ((mim2) 3 /(mi + m^)) 1 / 5 , 
and symmetric mass ratio r\ = mim?,/ (mi + m^) 2 ■ In the 
BBH search, we require coincidence in time, and the two phe- 
nomenological parameters and i/j^, which correspond to 
first approximation to M. c and r\ parameters, respectively (see 
i30ll46ll '). After the first inspiral filtering step, which does not 
use any computationally expensive vetoing methods such as a 
X 2 veto JUt], we apply coincidence windows with parameters 
that are summarized in TableHHI Then, in the PBH binary and 
BNS searches, we employ an hierarchical pipeline, in which 
coincident triggers are re-filtered, and the \ 2 veto is calcu- 
lated. Finally, trigger selection and coincidence requirements 
are re-applied. In the BBH search, no x 2 test is used because 
the waveforms have very few cycles in the LIGO detector fre- 
quency band. The coincident triggers from the first filtering 
step constitute the output of the BBH search. The coincident 
triggers from the second filtering step constitute the output of 
the PBH binary and BNS searches. 



TABLE III: Summary of the S3 and S4 coincidence windows. The 
second column gives the time coincidence windows; we also need to 
account for the maximum light travel time between detectors (10 ms 
between the LI and H1/H2 detectors). The third column gives the 
chirp mass (PBH and BNS searches), and ipo coincidence windows 
(BBH search). In the S4 BBH case, Atpo corresponds to about 1/15 
of the ipo range used in the template bank. The r\ (PBH and BNS 
searches) and ipz (BBH search) parameters (last column) are not 
measured precisely enough to be used in coincidence checks, except 
in the S4 BBH search. 





AT (ms) 


AM C (M ) 


A-q 


S3/S4 PBH 


4x2 


0.002 x 2 




S3/S4 BNS 


5x2 


0.01 x 2 






AT (ms) 


Aip 


Aip 3 


S3 BBH 


25 x 2 


40000 x 2 




S4BBH 


15 x 2 


18000 x 2 


800 x 2 



In the PBH binary and BNS searches, the \ 2 test provides a 
measure of the quality-of-fit of the signal to the template. We 
can define an effective SNR, p e g, that combines p and the \ 2 
value, calculated for the same filter, by 



2 

PeS 



f^J] d + 4- 

\2p~2 J \ 250 



(3) 



empirically, as justified in ||32Tl . We expect p e g ~ p for true 
signals with relatively low SNR, and low effective SNR for 
noise transients. Finally, we assign to each coincident trigger 
a combined SNR, p c , defined by 



(Pc)bNS,PBH 



N 



(4) 



where p is the number of bins used in the \ 2 test; the specific 
value of p = 16 and the parameter 250 in Eq. (01 are chosen 



where p c ff,i is the effective SNR of the trigger i th detector (HI, 
H2 or LI).' 

In the BBH search, no x 2 test is calculated. Therefore ef- 
fective SNR cannot be used. Furthermore, the combined SNR 
defined in Eq. (01 does not represent a constant background 
trigger statistic. Instead, we combine the SNRs from coinci- 
dent triggers using a bitten-L statistic similar to the method 
used in S2 BBH search IT30ll . as justified in l32tl . 

Finally, for each type of search, the coincident triggers are 
clustered within a 10 s window (BNS and BBH searches) or 
22 s window (PBH binary search), distinct from the clustering 
mentioned in Sec. Ill CI The final coincident triggers constitute 
the output of the pipeline — the in-time coincident triggers. 



III. BACKGROUND AND LOUDEST CANDIDATES 
A. Background 

To identify gravitational-wave event candidates, we need to 
estimate the probability of in-time coincident triggers arising 
from accidental coincidence of noise triggers, which consti- 
tute our background, by comparing the combined SNR of in- 
time coincident triggers with the expected background (with 
same or higher combined SNR). In each search, we estimate 
the background by repeating the analysis with the triggers 
from each detector shifted in time relative to each other. In the 
three searches, we used 50 time-shifts forward and the same 
number backward for the background estimation, taking these 
as 100 experimental trials with no true signals to be expected 
in the coincident data set. Triggers from HI were not time- 
shifted, triggers from H2 were shifted by increments of 10 s, 
and triggers from LI by 5 s. 

The time-shifted triggers are also used to explore the 
differences between noise and signal events in our multi- 
dimensional parameter space. This comparison is performed 
by adding simulated signals to the real data, analyzing them 
with the same pipeline, and determining the efficiency for de- 
tection of injected signals above threshold. This procedure 
allows us to tune all aspects of the pipeline on representative 
data without biasing our upper limits. The general philoso- 
phy behind this tuning process is not to perform aggressive 
cuts on the data, but rather to perform loose cuts and assess 
our confidence in a candidate by comparing where it lies in 
the multi-dimensional parameter space of the search with re- 
spect to our expectations from background. The details of this 
tuning process are described in detail in a companion paper 
rf32"ll . A representative scatter plot of the time-shifted triggers 
and detected simulated injections is shown in Fig.|2](S4 BNS 
case). This plot also shows how the effective SNR statistic, 
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• H1-L1 accidental events 
H1-H2 accidental events 

■ H1-L1 detected injections 
» H1-H2 detected injections 

• H1-H2-L1 detected injections 
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FIG. 2: Accidental events and detected simulated injections. This 
plot shows the distribution of effective SNR, p e ff, as defined in 
Eq. {3}, for time-shifted coincident triggers and detected simulated 
injections (typical S4 BNS result). Some of the injections are de- 
tected in all three detectors but no background triggers are found in 
triple-coincidence in any of the 100 time shifts performed. The Hl- 
Ll and H1-H2 time-shifted coincidence triggers have low effective 
SNR (left-bottom corner). 



which was used in the PBH binary and BNS searches, sepa- 
rates background triggers from simulated signals (with SNR 
as low as 8). 



B. Loudest candidates 

All searches had coincident triggers surviving at the end of 
the pipeline. In order to identify a gravitational wave event, 
we first compare the number of in-time coincident triggers 
with the background estimate as a function of p c . In S4, in- 
time coincident triggers are consistent with the background 
estimate in the three searches (see Fig. [3j. Similar results were 
obtained in S3 PBH binary and S3 BNS searches. However, in 
the S3 BBH search (not shown), one event clearly lies above 
expectation (in section UlIB 31 we explain why this candidate 
is not a plausible gravitational wave detection). The crite- 
rion we used to identify detection candidates which exceed 
expectation is to associate them with a probability P-q (p) that 
all background events have a combined SNR smaller than p. 
Pb(p) is calculated as the fraction of the 100 time-shifted ex- 
periments in which all triggers have smaller combined SNR 
than p. A candidate with a large Pb is considered a plausible 
gravitational wave event. If this is the case and/or a candidate 
lies above expectation we carefully scrutinize the data in the 
gravitational-wave channel and in auxiliary channels for pos- 
sible instrumental noise that could produce an unusually loud 
false trigger. We also investigate the astrophysical likelihood 
of the templates that best match the candidate in the different 
detectors (e.g., the ratio of effective distances obtained in dif- 
ferent observatories). In addition, irrespective of the outcome 
of the comparison between in-time and time-shifted coinci- 



dences, in-time coincident triggers with the highest p c values 
are also followed up. 

The loudest coincident triggers found in each of the 
searches are listed in Table [TVl Below, we briefly describe 
the reason(s) why we rejected the loudest candidates found 
in the three searches performed on the S4 run. These loudest 
events are used for the upper limit calculation (Sec. IIVI ). We 
also describe the loudest event found in the S3 BBH search 
mentioned above. 



1. Primordial black hole binaries 

There were no PBH binary candidates found in coincidence 
in all three detectors with SNR above the threshold of 6.5; 
nor were there accidental triple coincidences found in any of 
the 100 time-shifted runs. This means that had there been 
a triple-coincident candidate, there would be less than a 1% 
probability of it being a background event (Pb ^ 0.99). A 
cumulative histogram of the combined SNR of the loudest in- 
time coincident triggers in the S4 PBH search is shown in 
the leftmost plot of Fig. [3] The loudest S4 coincident trig- 
ger, with p c ff = 9.8, was found in coincidence in HI and LI. 
We observed equally loud or louder events in 58% of the 100 
time-shifted coincidence experiments. We found that this trig- 
ger was produced by a strong seismic transient at Livingston, 
causing a much higher SNR in the LI trigger than in the HI 
trigger; we found many background triggers and some missed 
simulated injections around the time of this event. As shown 
in Table IIVI the candidate also has significantly different ef- 
fective distances in HI (7.4 Mpc) and LI (0.07 kpc), because 
of the much larger SNR in LI: although not impossible, such 
high ratios of effective distances are highly unlikely. Tighter 
signal-based vetoes under development will eliminate these 
triggers in future runs. 



2. Binary neutron stars 

Just as in the PBH binary search, no triple coincident candi- 
dates or time-shifted triple coincident candidates were found 
in the BNS search. In-time coincident triggers were found 
in pairs of detectors only. We show in Fig. [3] (middle) the 
comparison of the number of coincident triggers larger than a 
given p c with the expected background for S4. The loudest 
coincident trigger was an HI -LI coincidence, consistent with 
estimated background, with p c = 9.1 and a high probability 
of being a background trigger (See Table llVI ). 



3. Binary black holes 

Due to the absence of a x 2 waveform consistency test, the 
BBH search suffered higher background trigger rates than the 
PBH binary and BNS searches, and yielded candidate events 
found in triple coincidence, both in S3 and S4. All triple 
coincident triggers were consistent with background. Never- 
theless, all triple coincidences were investigated further, and 
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FIG. 3: Cumulative histograms of the combined SNR, p c , for in-time coincident candidates events (triangles) and estimated background from 
accidental coincidences (crosses and 1 standard-deviation ranges), for the S4 PBH binary (left), S4 BNS (middle) and S4 BBH (right) searches. 
In each search, the loudest candidate (found in non-playground time) corresponds to an accidental coincidence rate of about 1 during the entire 
S4 run. 



none was identified as a plausible gravitational wave inspiral 
signal. In the rest of this section, we detail the investigations 
of the loudest triggers in each science run. 

In S4, the loudest coincident trigger in non-playground data 
was found in HI and H2, but not in LI, which was in operation 
at that time. This candidate has a combined SNR of 22.3 and 
Pb = 42%. The search produced many triggers in both HI 
and H2 at this time, reflecting a transient in the data produced 
by sharp changes in ambient magnetic fields due to electric 
power supplies. The magnetic fields coupled to the suspended 
test masses through the magnets used for controlling their po- 
sition and alignment. The transients were identified in voltage 
monitors, and in magnetometers in different buildings. The 
transients were rare, and were identified only in retrospect, 
when following up the loudest candidates, so they were not 
used as data quality vetoes in this analysis. 

In the playground data set, there was a louder candidate 
which has a combined SNR of 26.6 and Pb = 77%. This 
candidate was recorded during a time with elevated dust lev- 
els (due to proximate human access to the optics enclosure), 
which increases the transient noise in the detectors. Therefore 
this candidate was not considered to be a plausible gravita- 
tional wave event. 

In S3, the loudest candidate was found in coincidence in HI 
and H2, but not in LI, which was not in operation at that time. 
This candidate has a combined SNR of 107, resulting from a 
SNR of 156 in HI and 37 in H2. It lies above all background 
triggers and therefore has less than one percent probability of 
being background. None of the auxiliary channels of the Han- 
ford observatory show suspicious behavior at this time. This 
event was a plausible candidate and warranted further investi- 
gations via various follow-ups to confirm or reject a detection. 

We re-analyzed the segment at the time of this candidate 
with physical template families. At the coincidence stage, 
very wide coincidence windows in time (±25 ms) and chirp 
mass (±4 Mq), were required to get a coincident trigger. 
Then, based on the parameters of this coincident trigger, we 
analyzed the HI and H2 data around the candidate time with 
the same template. We compared the HI and H2 SNR time 
series; a real signal would produce a peak with the same time 
of arrival in both instruments to good accuracy. As seen in 



Fig. |H the maxima of both SNR time series are offset by 
38 ms, which is much larger than expected from simulations 
of equivalent gravitational wave waveforms with similar SNR 
and masses. Therefore, we ruled out this candidate from our 
list of plausible candidates. 

In summary, examination of the most significant S3 and S4 
triggers did not identify any as likely to be a real gravitational 
wave. 



IV. UPPER LIMITS 

Given the absence of plausible events in any of the six 
searches described above, we set upper limits on the rate of 
compact binary coalescence in the universe. We use only the 
results from the more sensitive S4 data and use only non- 
playground data in order to avoid biasing our upper limits 
through our tuning procedure. The u ppe r limit calculations are 
based on the loudest event statistic 1491 15011 . which uses both 
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FIG. 4: Time offset between the HI and H2 SNR time series, using 
the same template. Around the loudest candidate found in the S3 
BBH search (left panel), the maximum of the H2 SNR time series 
is offset by 38 ms with respect to the maximum of the HI SNR time 
series, which is placed at zero time in this plot. In contrast, simulated 
injections of equivalent gravitational wave waveforms with the same 
SNR and masses give a time-offset distribution centered around zero 
with a standard deviation about 6.5 ms (right panel). 
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TABLE IV: Characteristics of the loudest in-time coincident events found in the entire S4 data sets. Follow-up analysis of each of these events, 
described in section ITlIB I led us to rule them out as potential gravitational wave detections. Each loudest event was used in the final upper 
limit calculations. The first column shows the search considered. The second column gives the type of coincidence. The third column gives the 
combined SNR p c . The fourth column contains the parameters of the templates that produced the loudest triggers associated with this event. 
In the BNS and PBH binary searches, we provide the mass pairs mi , mi that satisfy coincidence conditions for chirp mass and symmetric 
mass ratio. The two masses can be significantly different because the coincidence condition on r\ is loose. In the BBH search, we provide the 
values of ipo and ^3. The fifth column is the effective distance in each detector which is provided for the BNS and PBH search only. The last 
column is the probability that all background events have a combined SNR less than p c . 





Coincidence 


Pc 


m 1 ,m 2 (M Q ) 


D off (Mpc) 


Pb(Pc) 


PBH 


(Hl-Ll) 


9.8 


(0.6,0.6) (HI), (0.9,0.4) (LI) 


7.4 (HI), 0.07 (LI) 


0.58 


BNS 


(Hl-Ll) 


9.1 


(1.6,0.9) (HI), (1.2,1.2) (LI) 


15 (HI), 14 (LI) 


0.15 


Vo(Hz 5 / 3 ),^ 3 (Hz 2 / 3 ) 


BBH 


(H1-H2) 


22.3 


(29000, -1800) (HI) 




0.42 


BBH 


(H1-H2) (playground time) 


26.6 


(153000, -2400) (HI) 




0.77 



100 




_ 10 20 30 40 50 

Effective Distance Livingston (Mpc) 

FIG. 5: Detection efficiency versus effective distance for the different 
searches (S4 run). The BBH and BNS efficiencies are similar, mainly 
because the loudest candidate in the BBH search is twice as loud as 
in the BNS search (See Table [TvTl. 



the detection efficiency at the combined SNR of the loudest 
event and the associated background probability. 

The Bayesian upper limit at a confidence level a, assuming 
a uniform prior on the rate R, is given by [f30h 



A 



1 + A 



i? T Cl (/O c ,max) 



(5) 

where Cl{p c . max) is the cumulative blue-light luminosity we 
are sensitive to at a given value of combined SNR p Ciinax , T 
is the observation time, and A is a measure of the likelihood 
that the loudest event is due to the foreground, and given by 



A 



\Pc,max) 



\Pc,max) 



pB\Pc,max) 



(6) 



where the derivatives are with respect to p c . As mentioned in 
Sec. Hill Pb (p) is the probability that all background events 
have a combined SNR less than p (shown in Table [IV] for 



the loudest candidates in each search). In the case where the 
loudest event candidate is most likely due to the background, 
A — > and the upper limit becomes 



#90% = 



2.3 



T \Pc,max) 



(7) 



In the limit of zero background, i.e. the event is definitely 
foreground, A — > 00 and the numerator in Eq. (0) becomes 
3.9. The observation time T is taken from Table HJ where we 
use the analyzed time not in the playground. 

The cumulative luminosity function Cl (p c ) can be obtained 
as follows. We use simulated injections to evaluate the effi- 
ciency £ for observing an event with combined SNR greater 
than p c , as a function of the binary inspiral chirp mass Ai c 
and effective distance D g. We then integrate £ times the 
predicted source luminosity L(D e s, A4 C ) as a function of ef- 
fective distance and mass. The detection efficiency is differ- 
ent for binary systems of different masses at the same effec- 
tive distance. Since we use a broad range of masses in each 
search, we should integrate the efficiency as a function of dis- 
tance and chirp mass. For low mass systems where the coa- 
lescence occurs outside the most sensitive region of the LIGO 
frequency band, the distance at which the efficiency is 50% 

5 /6 

is expected to grow with chirp mass: Z? o ff,50% M c (e.g., 
|47|P . We can define a "chirp distance" for some fiducial chirp 
mass M. c ,o as D c = D e g(A4 c ^ /A4 c ) 5 ^ ( ', and then measure 
the efficiency as a function of D c rather than D e g. This effi- 
ciency function is now independent of chirp mass, and the in- 
tegration can be performed with respect to the chirp distance 
only: Cl = J dD c £{D C ) L{D C ). We use a model based on 
101 for the distribution of blue luminosity in distance to cal- 
culate L(D C ) for a given mass distribution (e.g., uniform or 
Gaussian distribution). Since a system will have in general 
slightly different orientations with respect to the two LIGO 
observatories, they will also have slightly different effective 
distances. The efficiency for detection is thus a function of 
both distances, and the integration needed is two-dimensional: 
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Cl(p)= / / S(D c , H ,D c , L ,p)L(D c ,H,dD c>L )dD CiL dD c ,H. (8) 
Jo Jo 



The detection efficiency as a function of the effective distance 
for each observatory is shown in Fig. [5] This efficiency is 
computed using a Gaussian mass distribution, with a mean 
of M c ,o — 0.7Mq for the PBH binaries (mi = = 
0.75M o ), M c , ^ 1.2M Q for the BNS (m x = m 2 = 
1.4M Q ), M c ,o '~ 4.4M for the BBH (mi = m 2 = 5M Q ) 
and a 1M© standard deviation. These efficiencies are mea- 
sured with simulated injected signals, using the same pipeline 
we used to search for signals; the efficiency is the ratio of the 
number of injections detected with SNR above p c . max to the 
total number injected. We show in Fig.Q]the cumulative lumi- 
nosity as a function of effective distance in each observatory. 
It can be seen that the sharp drop in efficiency in Fig. [5] hap- 
pens at approximately the calculated horizon distance shown 
in Fig.[U 

The upper limit calculation takes into account the possible 
errors which arise in a search for PBH binaries and BNS, and 
are described in some detail in lIBlTl . We follow the analy- 
sis presented there to calculate the errors for the above result. 
The most significant effects are due to the possible calibra- 
tion inaccuracies of the detectors, (which are estimated by us- 
ing hardware injections), the finite number of Monte Carlo 
injections performed, and the mismatch between our search 
templates and the actual waveform. We must also evaluate 
the systematic errors associated with the astrophysical model 
of potential sources within the galaxy described in 0]. We 
obtain upper limits on the rate after marginalization over the 
estimated errors, as described in l5lll . 

In previous result papers (e.g., I27I0 . we used the Milky 
Way Equivalent Galaxy (MWEG) unit which is approximately 
1.7Lio, where Lio is 10 10 times the blue solar luminosity. 
In this paper, the merger rate estimates are normalized to 
galactic-scale blue luminosities corrected for absorption with 
the underlying assumption that merger rates follow the mas- 
sive star formation rate and the associated blue light emission. 
This assumption is well justified when the galaxies reached 
by the detector are dominated by spiral galaxies with ongoing 
star formation like the Milky Way. 

Assuming Gaussian mass distributions, as specified above, 
we obtain upper limits of 7Z 90 % = 4.9 yr -1 Lio" 1 for PBH 
binary, lZgo% = 1.2 yr -1 Lio" 1 for BNS, and lZgo% = 
0.5 yr -1 Lio 1 for BBH. We also calculated the upper lim- 
its as a function of total mass of the binary, from 0.7 Mq to 
80 Mq. These upper limits are summarized in Fig. [6] 

For comparison, we review the limits on the compact bi- 
nary coalescence rates from previous searches. The best pre- 
vious limits were obtained by TAMA300 using 2705 hours 
of data taken during the years 2000-2004; their result was 
an upper limit on the rate of binary coalescences in our 
Galaxy of 20yr~ 1 MWEG" 1 (12 yr~ x Lio" 1 ) JH- Previ- 
ous 90% limits from LIGO searches were 47 yr" 1 MWEG -1 
(28yr- 1 Li _1 ) in the mass range [1 - 3]M Q JH], and 
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38 yr" 1 MWEG" 1 (22 yr" 1 L^" 1 ) in the mass range [3 - 
20] M Q ll30ll (the numbers in brackets are in units per year per 
£10). 



V. CONCLUSION 

We searched for gravitational waves emitted by coalesc- 
ing compact binaries in the data from the third and fourth 
LIGO science runs. The search encompassed binary systems 
comprised of primordial black holes, neutron stars, and black 
holes. The search techniques applied to these data represent 
significant improvements over those applied to data from the 
second LIGO science run 11281 |29j, 13011 due to various signal 
consistency tests which have significantly reduced the back- 
ground rates at both single-detector and coincidence levels. 
Simulated injections with SNR as low as 8 are detectable, ex- 
tending the range of detection. In addition, the stationarity 
and sensitivity of the data from the S3 and S4 runs were sig- 
nificantly better than in S2. In the 788 hours of S3 data and 
576 hours of S4 data, the search resulted in no plausible grav- 
itational wave inspiral events. 

In the absence of detection, we calculated upper limits on 
compact binary coalescence rates. In the PBH binary and 
BNS searches, the upper limits are close to values estimated 
using only the sensitivity of the detectors and the amount of 
data searched. Conversely, in the BBH search, the short du- 
ration of the in-band signal waveforms and the absence of 
X 2 veto resulted in a significantly higher rate of background 
events, both at the single-detector level and in coincidence. 
Consequently, we obtained a reduced detection efficiency at 
the combined SNR of the loudest events and therefore a worse 
upper limit than we would have obtained using more effec- 
tive background suppression, which is under development. 
The upper limits, based on our simulations and the loudest 
event candidates, are 7^go% = 4.9, 1.2, and0.5yr _1 Lio" 1 
for PBH binaries, BNS, and BBH, respectively. These upper 
limits are still far away from the theoretical predictions (see 
Sec.[I]i. For instance, the current estimate of BNS inspiral rate 
is 10-170 x 10 -6 yr _1 ^ro 1 - 

We are currently applying these analysis methods (some- 
what improved) to data from LIGO's fifth science run (S5). In 
S5, all three detectors have achieved their design sensitivity 
and one year of coincident data are being collected. We also 
plan to use physical template families in the BBH search so 
as reduce the background and increase our confidence in de- 
tection. In the absence of detection in S5 and future science 
runs, the upper limits derived from the techniques used in this 
analysis are expected to be several orders of magnitude lower 
than those reported here. 
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FIG. 6: Upper limits on the binary inspiral coalescence rate per year and per Lio as a function of total mass of the binary, for PBH binaries 
(left), BNS (middle), and BBH (right) searches. The darker area shows the excluded region after accounting for marginalization over estimated 
systematic errors. The lighter area shows the additional excluded region if systematic errors are ignored. In the PBH binary and BNS searches, 
upper limits decrease with increasing total mass, because more-distant sources can be detected. In the BBH search, upper limits decrease down 
to about 30 solar mass and then grow where signals become shorter; this feature can be seen in the expected horizon distance as well (See 

Fig.CD- 
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